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摘 要 : 可 诱导 非 对 称 压 降 的 喷射 泵 被 引入 热 声 发 动机 的 环 路 结构 中 ， 可 以 用 于 抑制 Gedeon 流 ， 进 而 提高 热 声 发 动机 热 效 
率 。 针 对 Swift 对 喷射 泵 提供 时 均 非 对 称 压 降 的 理论 推导 结论 与 Idelchik 等 人 报道 的 实验 结果 之 间 存 在 的 矛盾 ， 本 文 利用 计 
算 流体 力学 软件 ， 研 究 了 喷射 泵 的 厚 径 比 Vs 对 时 均 非 对 称 压 降 的 影响 情况 。 结 果 表明 ， 改 变 厚 径 比 Va。 会 引起 时 均 非 对 称 
压 降 方向 的 变化 ,在 喷射 泵 设计 中 ,除了 应 该 考虑 喷射 泵 流 道 两 端面 积 不 相同 所 带 来 的 影响 ， 还 必须 合理 选择 其 厚 径 比 1/4。。 
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Analysis of Time-averaged Pressure Drop lnduced by Jet Pump 
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Abstract: Gedeon streaming, which may considerably deteriorate the thermal efficiency of a traveling-wave 
thermoacoustic engine, can be suppressed by inducing a jet pump into the loop configuration. The contradiction of 
the time-averaged pressure drop between Swift’s theoretical derivation and Idelchik’s experimental results, 
motivated us to study the effect of the jet pump’s thickness-to-diameter ratio Was on the time-averaged asymmetric 
pressure drop by means of computational fluid dynamics. The results indicate the variation in 
thickness-to-diameter ratio l/d. can change the direction of the time-averaged pressure drop induced by the jet 
pump, consequently, the designing of a jet pump needs to focus not only on the asymmetry of opening areas of 
both ends, but also on the thickness-to-diameter ratio /de. 
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0 引言 但 是 ， 反 馈 环 路 的 引入 会 在 热 声 发 动机 中 形成 

热 声 发 动机 由 于 具有 无 机 械 运动 部 件 、 结 构 简 闭合 环 路 , 这 就 可 能 导致 Gedeon 流 的 发 生 馈 。 所 谓 
证、 运 生 避 人 、 帮 全 长 和 特点 而 受到 人 们 的 广 污 关 。 。 Gedeon 流 是 沿 环 路 循环 流动 的 时 均 流 , 或 者 称 为 声 
道口， 理论 上 的 效率 高 于 驻 波 热 声 发 动机 ， 使 其 具 冷 端 换 热 器 而 不 参与 热 声 转换 过 程 ， 是 一 种 严重 的 
有 较 大 的 应 用 潜力 和 发 展 前 景 。Swif 等 人 研制 的 具 。。 朴 损 失 机 制 ， 其 成 为 影响 热 声 斯 竺 林 发 动机 热效率 


馈 环 路 的 执 声 斯 性 林 发 动机 D3] 旦 行 小 执 声 居 的 一 个 重要 因素 。 

有 反馈 环 路 的 热 声 斯 特 林 发 动机 231, 是 行 波 热 声 发 0 
动机 的 典型 代表 , 热效率 可 达 0.3, 引发 了 近年 来 关 Swift 等 人 采用 喷射 汞 ”对 Gedeon 流 进行 抑 
于 行 波 热 声 发 动机 的 研究 热潮 7。 制 。 利 用 两 端 流通 面积 不 相等 的 锥 形 槽 结构 ， 使 执 


声 发 动机 中 交 变 运动 的 流体 在 前 后 两 个 半 周 期 内 流 
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经 喷射 泵 时 ， 产 生 时 均 非 对 称 压 降 ， 进 而 起 到 抑制 
环 路 中 Gedeon 流 的 作用 。Swift 等 人 推导 的 相关 计 
算 公式 中 时 均 非 对 称 压 降 方 向 为 从 喷射 泵 的 大 截面 
到 小 截面 的 方向 ， 表 明 流 体 从 喷射 泵 的 大 截面 到 小 
截面 流动 的 压 降 大 于 反 向 流动 的 压 降 。 

Idelchik 等 人 外 研究 了 在 稳定 流动 的 情况 下 流 
体 分 别 以 相反 方向 流 经 非 对 称 孔 板 时 产生 的 压 降 情 
况 。 该 非 对 称 孔 板 与 喷射 泵 具有 相似 结构 ， 区 别 在 
于 Idelchik 等 人 报道 的 实验 研究 的 孔 板 为 圆 孔 而 非 
和 矩形 窗 槽 ， 且 孔 板 的 厚度 与 小 孔 的 直径 之 比 为 
1/4=0~0.015。 实 验 结果 表明 ， 流 体 从 小 截面 到 大 截 
面 流动 的 压 降 大 于 反 向 流动 时 的 压 降 。 

Swift 在 研究 喷射 泵 工作 机 理 时 ， 只 考虑 了 罕 
槽 流 道 两 端的 流通 面积 和 边沿 形状 不 同 这 两 个 因素 
所 造成 的 非 对 称 压 降 , 而 忽略 了 喷射 泵 厚度 的 影响 
这 可 能 是 导致 其 推导 结论 与 Idelchik 等 人 报道 的 实 
验 结 果 相 矛盾 的 原因 。 

基于 上 述 分 析 ， 本 文 参考 Swift 喷射 泵 模型 ， 
以 窄 权 型 喷射 泵 作为 研究 对 象 ,利用 数值 模拟 方法 ， 
研究 了 不 同 流 量 下 厚 径 比 Vd。 对 非 对 称 压 降 的 影 
响 ， 以 获得 适宜 的 二 阶 时 均 流 阻 且 保 持 较 小 的 一 阶 
交 变 流 阻 ， 达 到 高 效 低 损 抑制 Gedeon 流 的 目的 。 
1 计算 方法 的 验证 

参考 Swift 等 人 P3] 的 研究 思路 ， 将 交 变 流动 处 
理 为 沿 流 道 轴 向 的 两 个 相反 方向 的 稳定 流动 。 研 究 
在 稳定 流动 条 件 下 ， 喷 射 泵 厚 径 比 Wd 对 时 均 非 对 
称 压 降 的 影响 。 采 用 三 维 定常 不 可 压缩 雷诺 时 均 
N-S 方程 。 由 于 Realizable k-e 模型 能 描述 包含 有 射 
流 和 混合 流 的 自由 流动 、 管 道内 流动 、 边 界 层 流动 
等 不 同类 型 的 流动 09, 因此 在 模拟 中 选择 Realizable 
k-s (RKE) 模型 。 考 虑 到 在 喷射 泵 的 入 口 和 出 口 存 
在 较 大 的 压力 梯度 , 压力 的 插值 方法 选用 PRESTO。 
壁面 设 为 绝热 边界 条 件 , 流体 工 质 为 He。 连续 性 方 
程 、 动量 方程 、 能 量 方程 均 采 用 二 阶 迎 风 离 散 格 式 ， 
计算 时 认为 流体 在 管道 进出 口 的 流量 差 为 0.1% 时 
即 达到 收敛 .为 了 准确 模拟 喷射 泵 周围 的 涡流 信息 ， 
对 喷射 泵 以 及 邻近 直 管 段 进行 了 加 密 处 理 。 

为 了 验证 计算 方法 的 正确 性 ， 以 非 对 称 孔 板 为 
研究 对 象 外 ,对 流体 的 串 流 阻 力 系数 进行 了 计算 ， 
并 将 计算 结果 与 Idelchik 等 人 报道 的 实验 结果 进行 
比较 。 表征 流体 流 经 非 对 称 孔 板 产生 的 压力 损失 ， 
可 用 式 (1) 计算 : 


t=Ap/ 全 (1) 


2 
式 中 ,Ap 为 流体 流 经 非 对 称 孔 板 的 全 压 损 失 。 该 损 
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失 包 括 小 截面 处 的 截面 突变 引起 的 局 部 阻力 损失 、 
渐变 管 引起 的 摩擦 损失 、 浙 变 管 引 起 的 扩散 损失 以 
及 大 截面 处 的 截面 突变 引起 的 局 部 阻力 损失 四 部 
分 。p 为 流体 的 密度 ，* 为 管道 内 流速 。 

表 1 列 出 了 流体 从 非 对 称 孔 板 的 大 截面 到 小 截 
面 的 串 流 阻力 系数 的 计算 结果 与 实验 结果 的 比较 
情况 ， 其 中 管道 面积 与 非 对 称 孔 板 的 小 截面 的 面积 
比 4d/es 为 0.02，0.04，0.06。 管道 内 流速 vy 为 2 m/s。 
从 表 1 可 看 出 ， 计 算 结 果 与 实验 结果 的 误差 在 7% 
以 内 ， 可 见 该 计算 结果 能 够 合理 反映 实验 结果 。 

表 1 串 流 阻 力 系数 的 模拟 结果 与 实验 结果 比较 


Table 1 Comparison of simulation results and 


experimental data for the resistance coefficient 人 


面积 比 0.02 0.04 0.06 
实验 值 大 5157 1241.25 534 
模拟 值 4830 1163.07 504.61 

误差 -6.3% -6.29% -5.5% 


2 计算 参数 选择 

Swift 等 人 认为 时 均 非 对 称 压 降 仅 由 于 窜 槽 流 
道 两 端的 流通 面积 和 边沿 形状 的 不 同 引 起 。 本 文 为 
了 消除 以 上 两 个 因素 对 时 均 非 对 称 压 降 的 影响 ， 而 
究 在 喷射 泵 两 端的 流通 面积 不 变 以 及 边缘 均 不 倒 圆 
角 的 情况 下 , 厚 径 比 VC. 对 时 均 非 对 称 压 降 的 影响 ， 
并 在 此 基础 上 提出 合适 的 W4。.， 以 实现 在 提供 合适 
的 时 均 非 对 称 压 降 的 同时 使 压 损 保持 较 小 的 值 。 

采用 的 物理 模型 其 管道 内 径 为 60 mm， 喷 射 泵 
的 结构 参考 Swift 之 前 的 喷射 泵 模型 上 ]， 以 窗 槽 型 
的 喷射 泵 作为 研究 对 象 ， 其 小 截面 的 面积 为 as=1.8 
cm， 大 截面 的 面积 为 wu=3.0 cm2， 见 图 1。 计 算 采 
用 的 工作 压力 为 3 MPa， 温 度 为 300 K， 研 究 了 流 
速 为 2m/s，3 m/s，4 m/s 的 情况 下 ， 改 变 /as 对 时 
均 非 对 称 压 降 的 影响 ，4. 为 喷射 泵 的 小 截面 的 当量 
直径 。W4, 的 选择 : 0.2、0.5、1.0、1.5、2.0、3.0、 
4.5、5.5、6.5、7.0、7.5、9.0、10.0。 


1 和 矩形 截面 喷射 泵 简化 模型 
Fig. 1 Simplified model of jet pump with rectangular tapered 
channel 


3 计算 结果 

为 了 便于 分 析 ， 提 出 了 时 均 非 对 称 阻 力 系数 
AK、 时 均 总 阻力 系数 乒 和 有 效 系数 s 等 三 个 参数 用 
于 表征 喷射 泵 的 性 能 。 Ak 表征 流体 流 经 喷射 泵 产生 
的 时 均 非 对 称 压 降 ， 石 表征 流体 流 经 喷射 泵 时 产生 
的 时 均 总 压力 损失 ，s 表征 流体 流 经 喷射 泵 时 单位 
时 均 总 压力 损失 产生 的 时 均 非 对 称 压 降 。 它 们 的 定 
义 式 如 下 


7 7 
| Ap,dt—|,, ep @) 
下 - DY， 御 1 DY_ 1 有 
2 7/2 2 
| pdf Mp dt 
P+ 7/2 @) 
t 2 2 
| 全 dr 全 di 
0 2 7/2 2 
_Ak 
Fe @ 


式 中 ，7 为 周期 ， 下 标 +，-- 代 表 两 个 相反 的 流动 ，+ 
表示 流体 从 小 截面 流向 大 截面 ，- 表 示 流 体 从 大 截 
流向 小 截面 。 

图 2 为 流体 从 喷射 泵 的 小 截面 到 大 截面 方向 流 
动 的 串 流 阻 力 系 数 的 值 随 V4. 的 变化 。 图 3 为 流 
体 从 喷射 泵 的 大 截面 到 小 截面 方向 流动 的 串 流 阻 力 
系数 的 值 随 厚 径 比 V4, 的 变化 。 

从 图 2 可 以 看 出 ， 从 喷射 泵 的 小 截面 到 大 截面 
流动 的 串 流 阻力 系数 乒 随 厚 径 比 Wa. 的 增 大 而 逐渐 
减 小 ， 压 力 损失 减少 。 厚 径 比 Vd. 为 5.5~10.0 时 ， 
变化 趋势 逐渐 平缓 ， 随 着 厚 径 比 Wd. 进一步 增 大 ， 
压 降 有 略微 上 升 趋势 。 从 图 3 可 以 看 出 ， 喷 射 泵 的 
大 截面 到 小 截面 的 串 流 阻力 系数 左 的 变化 趋势 与 串 
流 阻力 系数 的 变化 趋势 相似 ， 但 大 的 变化 较 有 后 
的 变化 略 小 。 当 厚 径 比 W4, 大 于 3.0 时 ， 出 现 一 段 
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图 2 串 流 阻力 系数 丘 随 Wd. 的 变化 


Fig. 2 Variation in coefficient of resistance k; with ae 


图 3 是 


和 流 阻 力 系数 开 随 14. 的 变化 


Fig. 3 Variation in coefficient of resistance k_ with /de 


比较 平缓 的 曲线 ， 厚 径 比 V4. 的 增 大 对 串 流 阻力 系 


数 磊 的 影响 较 小 。 


图 4 和 图 5 显示 了 流速 为 2 m/s， 


厚 径 比 Wd 


分 别 为 0.2，2.0，5.5 时 流体 流 过 喷射 泵 的 速度 分 布 


情况 。 图 4 为 从 喷射 泵 的 小 
图 5 为 从 喷射 泵 的 大 截 
的 速度 分 布 图 。 

从 图 4 中 可 以 看 出 ， 对 了 
的 流动 ， 由 


度 分 布 图 ， 


大 截 
现 向 : 


质点 向 孔 口 加 速 ， 


三 面 到 大 截 


流动 的 速 


到 小 截面 流动 


六 从 喷射 
于 流体 受到 小 截面 的 约束 ， 会 出 


间 收 缩 的 情况 。 喷 射 泵 阻挡 流 
流体 流速 增加 ， 静 压力 降低 。 流 


束 进 入 小 截 


后 在 惯性 力作 | 


惯性 力 逐 渐 减 小 直至 与 扩张 力 平衡 ， 


Velocity 
Contour 1 


j 下 继续 收缩 ， 


泵 的 小 截面 向 


动 ， 迫 使 流体 


但 径 向 
这 时 流 束 的 横 


(b) 


.9 
:3 
.6 
.0 
] (c) 
图 4 不 同 Vd. 下 从 小 截面 到 大 截面 流动 的 速度 分 布 图 
(a) Vid.=02 (b) VA=2.0 (ce) Vd=5.5 


Fig. 4 Velocity distribution in flow from small opening to big 


(a) /d=0.2 


opening under various l/d。 
(b) Wd.=2.0 


(e) Vad=5.5 


与 Idelchik 等 人 报道 的 实验 结果 一 致 。 随 着 Vi。 的 


Co 增加 ， 流 体 从 喷射 泵 的 小 截面 流向 大 截面 的 压 降 迅 
38 速 降低 , 而 反 向 流动 的 压 降 减 小 相对 较 慢 ， Ak 随 之 
[34 减 小 ， 当 Vad, 进一步 增 大 ， 流 体 从 喷射 泵 的 小 截面 
:30 (a) 流向 大 截面 的 压 降 最 终 小 于 反 向 流动 压 降 , 此 时 Ak 
26 变 为 负数 .Swift 理论 模型 关于 喷射 泵 时 均 非 对 称 压 
E32 降 的 方向 预测 与 Ak 变 为 负数 时 的 情况 一 致 ， 而 与 
得 Ak 为 正 数 时 的 情况 相反 。 究 其 原因 , 主要 在 于 Swift 
16 加 等 人 在 进行 理论 分 析 时 ， 忽 略 了 两 个 突变 截面 局 部 
12 阻力 间 的 相互 影响 。 

150 
4 
2 
O 
[m s^-1] 


(9) 


一 5 不 同 Via。 下 从 大 截面 到 小 截面 流动 的 速度 分 布 
PP (a Vds=0.2  (b) Wd.=2.0 (©) Ud=5.5 
Fig. $ Velocity distribution in flow from big opening to small 
opening under various /al。 
(a) UVds=0.2 (b)Wd=2.0  (c) Vd=5.5 


| 


=。 ”截面 积 达 到 最 小 , 流速 增加 到 最 大 , 静 压 降 到 最 低 。 Ud. 
_ ”此 后 ， 流 束 扩 大 ， 速 度 降低 ， 直 至 流 束 充满 管道 。 6 时 均 非 对 称 阻力 系数 Ak 随 /ix 的 变化 
| 在 Wd. 较 小 的 情况 下 ， 由 于 截面 突 缩 带 来 的 流 Fig. 6 Variation in coefficient of time-averaged resistance Ak 
二 束 收 缩 有 一 定 影响 长 度 ， 在 计算 出 口 膨胀 的 压力 损 es 


SJ] ” 失 时 ， 必 须 考虑 在 边 壁 突变 的 附近 出 现 的 主流 与 边 
>” 壁 脱 离 的 现象 。Vie 为 5.5 时 ,实际 流 束 面积 接近 于 
一 出 口 壁面 面积 ， 主 流 与 边 壁 脱离 的 影响 减 小 。 这 时 
个 串 流 阻力 系数 大 的 波动 ， 主 要 由 喷射 泵 内 部 摩擦 损 
失 与 扩散 损失 随 厚 度 的 变化 所 引 
r= 从 图 5 中 可 以 看 出 ， 流 体 从 喷射 泵 的 大 截面 到 
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图 7 为 时 均 总 阻力 系数 的 变化 情况 ， 其 随 省 

yd 增 大 而 降低 ,之 后 变化 趋 组 ， 随 着 74, 进一步 增 

大 ， 总 压力 损失 咯 微 上 升 

图 8 为 有 效 系数 。 的 示意 图 ， 为 了 实现 高 效 低 

损 的 抑制 直流 ，& 的 值 应 尽 可 能 大 。 从 图 中 可 以 看 

区 小 由 和 增 大 呈现 先 减 小 后 增 

小 截面 流动 时 ， 由 于 鹤 面 为 渐 缩 ， 流 体质 点 受到 与 。。 天 生 区 仆人 
流动 方向 一 到 的 正 压 差 作用 ， 只 会 导致 加 速 ， 不 会 

出 现 主 流 与 边 壁 脱离 的 现象 。 500 

从 模拟 结果 可 以 看 出 ， 管 道内 流速 的 改变 对 
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流 阻力 系数 的 影响 较 小 。 本 文 将 交 变 流动 处 理 为 沿 
流 道 轴 向 的 两 个 相反 方向 的 稳定 流动 ， 假 定 前 后 半 本 
个 周期 的 流速 大 小 不 变 。 因 此 ， 定 义 式 (2) 和 (3 ) 本 
可 简化 为 350 
es 太一 大 (5) 300 
2 2 2 
站 250 
本 Ap +Ap_ _ k +k (0) 
2 2 2 1 
图 6 为 时 均 非 对 称 阻力 系数 Ak 随 Wd 的 变化 。 图 7 时 均 总 阻力 系数 乒 随 Vd. 的 变化 
从 图 中 可 以 看 出 ， 1 较 小 时 ， 流 体 从 喷 射 泵 的 小 Fig. 7 Variation in coefficient of overall resistance ki with az。 


截面 流向 大 截面 时 的 压 降 大 于 反 向 流动 的 压 降 ， 这 
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治 相 反方 向 流动 的 压 降 ， 时 均 非 对 称 压 降 的 方向 为 
从 喷射 泵 的 小 截面 到 大 截面 的 方向 , 其 结果 与 Swift 
等 人 理论 模型 的 推导 结果 相 矛 盾 。 随 着 Vae 的 增加 ， 
流体 从 喷射 泵 的 小 截面 流向 大 截面 的 压 降 迅速 降 
低 ， 最 终 小 于 反 向 流动 压 降 ， 时 均 非 对 称 压 降 的 方 
向 为 从 喷射 泵 的 大 截面 到 小 截面 的 方向 , 这 与 Swift 
理论 推导 的 结论 一 致 。 厚 径 比 Vd 改变 ， 使 时 均 非 
对 称 压 降 的 方向 发 生变 化 。 因 此 , 在 设计 喷射 泵 时 ， 
必须 综合 考虑 两 端面 面积 和 厚 径 比 Vae 的 影响 ， 否 
则 可 能 会 导致 喷射 泵 的 设计 不 合理 ， 不 仅 不 能 抑 于 
直流 ， 甚 至 会 在 一 定 程度 加 剧 环 路 中 的 直流 。 
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